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The ground and excited vibrational states for the three hydrogen isotopes 
on the Pd(lll) surface have been calculated. Notable features of these 
states are the high degree of anharmonicity, which is most prominently 
seen in the weak isotopic dependence of the parallel vibrational transition, 
and the narrow bandwidths of these states, which imply that atomic 
hydrogen is localized on a particular surface site on time scales of 100 
picoseconds or more. Experiments to resolve ambiguities concerning the 
present system are suggested. 
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Hydrogen atoms chemisorbed on a metal surface [1] display a variety of 
phenomena, ranging from localization on a particular binding site with 
thermally activated hopping or incoherent tunneling between the sites to 
extended band states, depending on the metal, the face, the tempera- 
ture, and the coverage. A hydrogenic-band model in which the hydro- 
gen wavefunction is delocalized parallel and localized perpendicular to 
the surface was first [2] proposed by Christmann, et a/., in 1979 for H 
on Ni(lll)[3]. Puska and co-workers followed up on this suggestion by 
calculating the states for hydrogen on a rigid Ni surface [4]. They find 
narrow width ground states and excited states with appreciable width 
(~50 meV), particularly for the (110) and (111) face. Since then there 
has been a number of experiments which reported evidence of band-like 
states, mainly along the "troughs" of surfaces such as the (211) face of 
body-centered cubic (bcc) metals[5] and the (110) face of face-centered 
cubic (fee) metals [6] or for the smoothly corrugated (111) face of fee met- 
als[7,8]. The helium scattering experiments of Hsu, et al. find a hydrogen 
phase with 6- fold symmetry [8]. The authors explain this phase in terms 
of hydrogen being quantum mechanically delocalized over two adjacent 
3-fold sites on the sub-picosecond time scale of the helium scattering 
event. These results are for a hydrogen coverage of half saturation and 
a temperature range of 140 to 300 K and the same results are seen for 
D as for H [9]. The time scales for hydrogen diffusion have been mea- 
sured for several metal surfaces (Ni(100)[10],W(110)[ll],Rh(lll)[12], and 
Pt(lll)[13]). At high temperatures, diffusion shows an Arrhenius depen- 
dence indicative of thermally activated hopping between localized sites, 
with activation energies ranging from 0.15 to 0.52 eV. At temperatures 
below 100 K, Gomer and coworkersfll] for H/W(110) and Zhu, et al.[U], 
for H/Ni(100) find a diffusion rate due to tunneling which is independent 
of temperature. The diffusion rates at the low temperatures are very 
small (10 -11 to 10~ 13 cm 2 /sec), indicating very little wavefunction over- 
lap between adjacent sites and bandwidths of 10~ 2 to 10 -11 meV [15]. 
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The small bandwidth for H on Ni(100) is consistent with the band cal- 
culations of Puska and co-workers [4] and high resolution electron-energy 
loss spectroscopy (HREELS) experiments[16]. 

Vibrational spectroscopy such as HREELS provides a direct exper- 
imental probe of the nature of the hydrogen states. HREELS gives di- 
rect evidence as to whether vibrational excitations correspond to small 
vibrations about an equilibrium position as is typically the case (e.g. 
H/Ni(100)) or correspond to transitions to broad bands (H/Cu(110)). An 
HREELS study of H/Pd(lll) at high coverages (a hydrogen exposure of 
2 Langmuirs) shows two narrow peaks: a parallel vibration at 774 cm -1 
and a perpendicular vibration at 998 cm -1 [17]. In this letter we report 
the vibrational states for the three isotopes of hydrogen on Pd(lll) in 
the single adatom limit. The low coverage limit presents the best condi- 
tions for delocalized states, since the H-H interaction is expected to be 
purely repulsive. Increasing coverage should lead to a higher probabil- 
ity of localized states. This supposition is borne out by experiment, in 
which evidence of quantum derealization disappears at high coverages 
[6,8]. Furthermore, although the widths of the excitations can vary with 
coverage, the position of the peak center does not have a strong coverage 
dependence. 

In the present computations, the lattice is kept rigid in the minimum 
energy configuration of the surface in the absence of any adsorbates. 
The magnitude of adsorbate induced reconstructions is examed below. 
The clean surface does not reconstruct, the top surface only contracts 
by 0.07 A relative to the bulk lattice spacing. The (111) surface is a 
hexagonal surface with 2 distinct sites: a tetrahedral site directly above 
a subsurface atom and an octahedral site, over a metal atom in the 
second subsurface layer. The two sites are treated as fully degenerate 
since the potential used predicts them to be energetically equivalent. 
The interaction potential is taken to be the Embedded Atom Method 
(EAM) potential of Daw and Baskes [18]. The EAM potential has a form 
borrowed from density functional theory and is determined by fitting to 
properties of the bulk system. Although this potential contains only bulk 
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data, it successfully reproduces the phase diagram for H/Pd(lll)[19] 
and predicts the heat of absorption of hydrogen to be 2.85 eV, which is 
the experimental value [20]. This predicted value is the classical binding 
energy relative to an isolated atom and is for hydrogen absorbed on the 
rigid surface. Lattice relaxations lower the energy an additional 0.077 
eV. The EAM barrier for surface diffusion is 175 meV. This barrier is not 
known experimentally but the barrier to bulk diffusion is 226 meV[21]. 
EAM predicts 152 meV for the bulk barrier, allowing for relaxation of 
the palladium atoms. 

To find the vibrational states, we use the tight-binding approxima- 
tion in which we calculate states localized on a particular binding site 
and take the appropriate linear combinations to construct states with the 
full 2-dimensional periodicity of the surface (see, e.g., Ref. [22]). The full 
periodic Hamiltonian of our system, 7i, is writen as a localized part, TIq, 
plus the corrections, AU. In the localized Hamiltonian, the hydrogen only 
interacts with the nearest three surface palladium atoms and the nearest 
three subsurface palladium atoms. These six interactions most nearly 
reproduce the full potential of the lattice in the region of the 3-fold site. 
Three eigenstates of Ho are calculated: a ground state with Ai symme- 
try and excited states with Ai and E symmetries. The states are found 
by solving the Schrodinger equation on a 3-dimensional mesh of spacing 
0.025 A and — for the the excited states — simultaneously projecting out 
the lower energy states [23]. The accuracy of the ground and E states 
can be checked using the Diffusion Monte Carlo (DMC) method, which 
solves the Schrodinger equation exactly [24]. The E state can be calcu- 
lated using the DMC method using the fixed-node approximation, since 
for the E-symmetry state the fixed node approximation is exact and the 
position of the nodal plane is easily deduced by symmetry. For the Ai 
state, the fixed-node approximation is not exact. The Bloch states, \I/n,k, 
with k the two-dimensional wavevector, are constructed by calculating 
the overlap between nearest neighbor surface sites only. The necessary 
overlap integrals can be easily done on the same mesh as is used to solve 
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for the localized states. The overlap integral between adjacent sites is 
less than 0.17, so the tight-binding approximation is reasonable. 

The band centers and widths for the three lowest energy states are 
shown on Table I. The energies are all relative to the classical binding 
energy (-2.85 eV). The states for all three isotopes are very localized, 
the only states with a bandwidth larger than 1 meV are the H excited 
states. The width for each state narrows with increasing mass roughly 
as e-^™[25]. 

In order to assist with spectroscopic assignments the k=0 H band- 
states are plotted along the three directions (Fig. 1). The left side of Fig. 
1 shows the states' amplitude as a function of one surface coordinate, the 
[121] direction, and the [111] direction perpendicular to the surface, at 
a constant value for the other surface coordinate equal to the center of 
the unit cell. The right side shows the same wavefunctions as a function 
of the two surface coordinates at a constant height of lA above the sur- 
face plane. The localized nature of the ground state is apparent. The 
excited state A} is much broader and can be identified as predominantly 
a perpendicular vibrational excitation, since the nodal surface is nearly a 
plane parallel to the surface. The to A{ transition will be dipole ac- 
tive in HREELS, the A^ to E transition can be seen only though impact 
scattering [26]. 

Table I also shows the energy differences between the ground and 
excited states for the parallel mode (the A° to E transition) and the 
perpendicular mode (A^ to A}) for the three isotopes. Shown is the en- 
ergy differences between the centers of the bands, although for the H 
isotope the excited states are broad enough that wide loss peaks might 
be detectable experimentally. The ordering of the vibrational peaks is 
the opposite of Conrad, et a/., who assign cu||=774 cm -1 and uj±=998 
cm -1 based on their HREELS experiments at high coverages [17]. The 
ordering of the parallel and perpendicular modes does vary for hydrogen 
on various metals and is a sensitive gauge of the metal-hydrogen inter- 
action[l]. For examples with three-fold or quasi-three-fold coordination, 
H/Ru(0001) [27] and H/Ni(lll) [28] have u± > u\\ and H/Pt(lll) [29] 
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and H/Rh(110)[30] have u\\ > uj±, although the H/Pt(lll) assigment is 
challenged by Feibelman and Hamann [31]. If the experimental assign- 
ments for H/Pd(lll) are correct, the EAM potential used here appears 
to misrepresent the relative stiffness of the two vibrations. 

The isotope dependence of the computed frequencies is unusual and 
demonstrates the importance of anharmonic effects, particularly for the 
parallel frequency. The ratio of uj\\ (H)/u;|| (D) is only 1.19, different from 
the a/2 dependence of harmonic vibrations. The ratio of u±(ii)/u>±(D) 
is 1.32, closer to \/2. The anharmonic nature of the vibrational states 
can be seen by comparing the exact energies with harmonic estimates 
(Fig. 2). The band state energies (the longer horizontal line on the 
left) is shown connected (by the dotted line) to the harmonic estimate of 
the energy of that state (the smaller horizontal line on the right). The 
harmonic estimate is based on the curvature of the potential minimum 
of the absorbed state. Also visible in Fig. 2 are the bandwidths of the 
A{ and E states. All the states are lower in energy than the harmonic 
estimate. Most notable is the E state of H, which is lower by 90 meV 
or about 20%. That the vibrational frequencies are sensitive to features 
of the potential other than the curvature of the well bottom should be 
a warning to those who might attempt to construct potentials by fitting 
harmonic estimates to experimental values of the frequencies. 

Most experimental spectra have a isotope dependence closer to y/2, 
but the weak dependence of the parallel mode has precedence: HREELS 
experiments by Mate and Somorjai for H/Rh(lll) find only a very small 
isotopic shift of about 10 cm -1 for a loss feature believed to be a transition 
to a delocalized band with E symmetry [7]. HREELS experiments for 
H and D on Ni(100)[32] and Rh(100)[33] have frequency ratios of the 
perpendicular mode close to y/2, although overtones in the H/Rh(100) 
spectra show noticeable anharmonicities. For H and D on Pd(100) [34], 
the uj±_ ratio is about 1.3 and for the Ni(110) and Rh(110) surfaces the 
ratio are 1.25[35] and 1.35 [30], respectively. Calculations for H/Rh(0001) 
find isotope ratios for several states all in the range of 1.41 to 1.46 and 
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harmonic estimates differ from the anharmonic calculations by, at most, 
8% [36]. 

The H/Pd(lll) bandwidths are much narrower and the band gaps 
larger than those calculated for H/Ni(lll)[4]. The ground state for that 
system has a width of 4 meV, the first excited state, which has E sym- 
metry, has a width of 40 meV, and the next excited state, with Ai sym- 
metry, has a width of 74 meV. The larger widths of the Ni states is due 
to the lower diffusion barrier for H on Ni(lll)[18] and the shorter lattice 
constant for Ni (3.52A for Ni, 3.89 A for Pd). 

From the localized tight-binding states, estimates of both the mag- 
nitude of hydrogen induced reconstructions and time scales for hydrogen 
tunneling between adjacent surface sites can be made. As mentioned 
previously, lattice distortions lower the classical binding energy by 77 
meV, as predicted by the EAM potential. This energy is a factor of ten 
larger than the largest bandwidth (Table I) and four orders of magnitude 
greater than the ground state width. Thus the small polaron condition is 
easily met: the energy gained from self-trapping the adsorbate by small 
lattice distortions is greater than the kinetic energy gained by forming ex- 
tended bands[37]. Small polaron theory provides a convenient formalism 
for describing the tunneling process from site to site, including the influ- 
ence of phonons. The small polaron Hamiltonian is found by expanding 
the potential around the rigid lattice potential, V r * s * d (r), to first order 
in the phonon coordinates, 



V s P(r, Q) = V r ^ d (r) + ^ VfQi 



(1) 



where 



V? = 



<9V(r,Q) 



<9Qi 



V(r, Q) is the EAM potential energy, Qi is the I th perfect lattice (no 
adsorbate) phonon coordinate and r a is the hydrogen position of the 
minimum of V ri5ld (r) at site "a". The small polaron Hamiltonian, H p , is 
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then the sum of the rigid lattice Hamiltonian, ri, a perfect lattice phonon 
part and the linear coupling between the two, 

Hp = n + J2 v? + \m Pd ^Q 2 ^ + Vi a Qi 

where mpd is the palladium mass and u)\ the I th phonon frequency. The 
H-phonon coupling will cause the phonon coordinates to be displaced by 
V L a /mpdujf. The energy gained by the phonon displacements is E\ oc = 
£(^)72mp^ 2 . 

In the temperature regime where tunneling is the dominant mecha- 
nism for diffusion (T< 150-100 K) only the ground state will be populated 
due to the large gap between the ground and excited states. The time 
scale, T a b, for tunneling from the ground state of site "a" (0 O ) to any of 
the three nearest neighbor states "b" {4>b) is given by 

r~ h X = 3J ab S ab (T)/h. 

and is the product of the Hamiltonian matrix element between adjacent 
tight-binding states, J a b = ((f>a\'Hp\(f>b) , and a temperature dependent 
phonon overlap term, S ao . The factor of three accounts for the number 
of nearest neighbor sites and the matrix element J ab is equal to 1/6 
the groundstate bandwidth [38]. Phonon coupling destroys some of the 
coherence of the bands, an effect which increases with temperature. At 
T=0, 

S ab (T = 0) = e -£. M^)((^"^)/™) 2 . 

Just considering the 24 phonon modes from the vibrations of the four 
atoms involved in the unit cell plus the four nearest neighbor surface 
atoms gives a E ioc of 46 meV. The localization energy is more than half 
the exact localization energy (77 meV) gained by relaxing the entire lat- 
tice, so the 24 phonon modes included here contain most of the impor- 
tant relaxation modes. The phonon displacements that give rise to this 
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localization are small. For the most strongly coupled mode (a symmetric 
stretch parallel to the surface of the three nearest Pd atoms) , the phonon 
displacement, Vf/mpawf, is only 0.08 A. 

The time scale for tunneling from Equation (3) is 2.6xl0 2 ps for H, 
3.8xl0 4 ps for D, and 7.2xl0 5 ps for T. The phonon overlap, S a b(T=0), 
is 0.46 and is only small modification on the tunneling matrix element. 
These time scales are all much longer than the sub-picosecond time scale 
needed to explain the Hsu, et a/., experiments within a quantum der- 
ealization mechanism [8]. The calculated rates are the result of a model 
potential and are by no means exact, but there are a few reasons to be- 
lieve that the tunneling time scales should be greater than 1 ps. Since 
both the diffusion barrier and the lattice constant for the Ni(lll) surface 
are less than for Pd(lll), the H/Pd(lll) bands should be narrower than 
4 meV, the calculated H/Ni(lll) groundstate bandwidth [4]. This band- 
width, when inserted into Equation (3), gives a time scale of .3 ps (with- 
out phonon corrections); a narrower bandwidth would give localization 
on the helium scattering time scale. In addition, since the isotope depen- 
dence is so strong, the experimental results, if due to large bandwidths, 
should be different for H and D. The evidence for derealization from 
the Hsu experiments is the symmetry change with decreasing coverage 
from 6-fold to 3-fold and the accompanying attenuation of the specular 
intensity. The delocalized states are deduced using a model for the he- 
lium scattering process in which the helium atoms scatter off a hydrogen 
density which in the plane of the surface is a constant value within an 
ellipsoid centered at the bridge site and a major semi-axis equal to 0.8A 
and minor semi-axis equal to 0.5A. Thus the hydrogen density has a 
constant value which spans both sides of the unit cell and the hydrogen 
is delocalized within the unit cell. This is a qualitatively different picture 
of the hydrogen wavefunction than the ground state shown in Fig. 1. 

There are factors other than the accuracy of the EAM potential 
which are relevant in a comparison of the present calculations with the 
experiments of Hsu, et al. The experiments were done at a coverage of a 
half so hydrogen-hydrogen interactions might be an issue, although most 
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evidence supports the conclusion that low coverage presents the best 
conditions for derealization (see above). Another issue is the presence of 
subsurface sites, which are of comparable energies with the surface sites 
and are thermally accessible. These sites, which are not visible to helium 
scattering, might influence the symmetry of the surface phase[19]. 

In summary, the vibrational states for H,D, and T on the Pd(lll) 
surface were reported. A notable feature of these states is their anhar- 
monicity, which is most prominent in the weak isotopic dependence of 
the parallel vibrational state (see Table I). The poorness of harmonic es- 
timates of the energies also demonstrates the importance of anharmonic 
effects. Harmonic estimates of the energy of the E state of H are off by 
90 meV. Harmonic estimates do better for the A{ state and get better 
with increasing mass (see Fig. 2). The states have a large zero point 
energy, so it is not surprising that anharmonic regions of the potential 
energy surface are important. The ground state energy for H and the 
excited state energies for all three isotopes (see Table I) are greater than 
the classical barrier between sites, 175 meV. That these states can have 
a high energy relative to a barrier in one coordinate and still remain 
largely localized on one side of this barrier is due to the multidimen- 
sional nature of these states. For a discussion of the localization and zero 
point/barrier energy ratio in the context of clusters see Ref. [37]. The 
bandwidths of the states are narrow and an analysis of the coupling of 
the hydrogen eigenstates to the phonons indicates that the hydrogen only 
induces small lattice distortions (with an amplitude < 0.1 A) and that 
the hydrogen is localized on a particular binding site for a time scale on 
the order of a hundred picoseconds or more. The dynamics of hydrogen 
as it diffuses across the surface and also into the subsurface will be the 
subject of future study [40]. 

That anharmonic effects are more important for the computed per- 
pendicular vibration (transitions to the E state) than the parallel vibra- 
tion (transitions to the A{ state) could lead to some atypical HREELS 
results. HREELS finds that the the parallel frequency, uju, is less than 
the perpendicular, uj±, for the H isotope[17]. Since the perpendicular 



10 



mode has a stronger isotopic dependence than the parallel, it is possible 
that the ordering of the frequencies could change for heavier hydrogen 
isotopes and u\\ could become higher than u±, although the present cal- 
culations find uj± < ui\\ for all isotopes. An HREELS study of H and 
D on Pd(lll) could resolve some of the issues discussed in this letter, 
including the ordering of the parallel and perpendicular vibrations and 
the isotopic dependence of these frequencies. Furthermore, the HREELS 
spectra could determine whether the excited states have bandwidth which 
is much larger than HREELS resolution (about 1 meV) or of the same 
magnitude as we are predicting here. 
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Table I. Band centers and band widths for the ground and two first 
excited vibrational states of the hydrogen isotopes on the Pd(lll) surface, 
in meV, and the vibrational frequencies corresponding to the the to 
E transition (u\\) and the A° to A \ transition (u>±), in cm -1 . 



H 




D 




T 




Center 


Width 


Center 


Width 


Center 


Width 


(meV) 


(meV) 


(meV) 


(meV) 


(meV) 


(meV) 


A? 236.3 


0.011 


169.2 


0.000076 


138.6 


0.000004 


A\ 345.9 


2.17 


252.5 


0.023 


208.2 


0.0006 


E 356.7 


7.7 


270.7 


0.208 


228.35 


0.0075 


Frequencies 


cuii/cm -1 971 




819 




724 




wi/cm" 1 884 




672 




561 





Figure Captions 

Figure 1. Contours of the amplitudes of the k=0 bandstates for H/Pd(lll) 
as a function of a surface coordinate and a coordinate perpendicular to 
the surface (left side) and as a function of the two surface coordinates 
(right side), for the E (top), A{ (middle), and A? states (bottom). Each 
rectangle spans an area of 2.75 A by 4.76 A and the coordinate in the 
[111] direction begins at .02 A below the surface plane. 

Figure 2. Bandstate energies for H, D, and T, comparing the tight- 
binding state energies (the longer horizontal lines) to the harmonic esti- 
mates (the shorter horizontal lines). 
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